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NATIONAL ADVISCET COMMITTEE FOB AEEONADTICS 


TECHNICAL NOTE 2295 


CHORDVn:SE AND COMERESSIBILITT CQREECTIONS 
TO SIENDEE-WING THEORY 
By Harvard Lomax and Loma Sluder 


SUMMARY 


Corrections to slender-wing theory are obtained by assuming a span- 
wise distribution of loading and determining the chordwise variation 
which satisfies the appropriate integral equation. Such integral equa- 
tions are set up in terms of the given vertical induced velocity on the 
center line or, depending on the type of wing plan form, its average 
value across the span at a given chord station. The chordwise .distribu^ 
tion is then obtained by solving these integral equations. Eesults are 
shown for flat— plate, rectangular , and triangular wings. 


INTRODUCTION 


The calciilation of loading on three-dimensional lifting surfaces is 
a fundamental problem in aerodynamic research. . The complexity of the 
problem has lead to the development , of certain simplified theories by 
means of which the loading on special types of plan forms can be esti- 
mated quickly. The amount of error which these estimates contain is of 
considerable interest, as are methods which will tend to correct such 
errors without undue labor. 

Slender-wing theory applies to one such sinplified body of analysis. 
There are two basic assumptions of this theory: one, the angle of attack 
is small enough so that the vortex sheet does not separate from the wing 
and the boimdary conditions for the wing can be projected onto a hori- 
zontal plane parallel to the direction of the free stream; and the other, 

. that either the chordwise gradient of velocity is small enough or the 
free— stream Mach number is close enoiigh to unity that the linearized 
partial differential equation which governs the fluid flow becomes 
Laplace’s equation in a plane transverse to the free— stream direction. 
Eeferences 1 through 5 axe examples of recent papers developing slender- 
wing theory. 
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An indication of the magnitude of the error which is introduced hy 
the use of such a theory comes, in the subsonic case, by observing that 
solutions so obtained violate the Kutta condition at the trailing edge. 
Proper inclusion of the chordwise and compressibility effects must 
result in solutions which satisfy the Kutta condition and make the load- 
ing, fall to zero at the trailing edge. It is the purpose of this report 
to study such modifications. 

The corrections due to the chordwise and compressibility effects 
are obtained in the following manner: First, an integral equation involv- 
ing elementary horseshoe vortices is set up; second, this integral equa- 
tion is solved under the assumption that the chordwise gradients are 
small or that the free— stream Mach niunber is \mity; and finally, the 
integral eqiiation is reinspected, this time with the spanwise lift dis- 
tribution fixed at the shape just obtained and with the chordwise varia- 
tion as the unknown eind the Mach number terms included. 

Results are presented and discussed both for triangular and rectan^ 
gular, flat— plate plan forms in both subsonic eind supersonic flow. 


LIST OF IMPORTANT SYMBOLS 


Ar 

b 


'm 


E 


velocity of sound in the free stream 
aspect ratio (f) 

reduced aspect ratio (Aj. = 3A) 

span of wing measxired normal to plane of symmetry 


root chord of wing 


/ lift 
V qS 


) 


lift coefficient 
pitching-moment coefficient 


f pitching moment 

qSCo 




) 


complete elliptic integral of the second kind with modulus k 

^l-q£2t2 


t A-k^t" 



E(t,k) incomplete elliptic integral of the second kind with argument t 

E(t,k) = 


and modulus k 



dt 


l-t2 
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F(t,k) incomplete elliptic integral of the first kind with argument t 


and modulus k F(t,k) = 


'o 7(l-t2)(i-^s2t2) 


complete elliptic integral of the first kind with modulus k 




for triangular wing, slope of leading edge relative to plane of 
symmetry 


free— stream Mach number 


static pressure 


(S) 


AP Pz - Pu 

q free— stream dynamic pressure ^ I PoVo^) 

s semispan of rectangiilar wing 

S area of wing 

u perturbation velocity component in the direction of the x axis 

Au Uu - uj 

Vq free— stream velocity 

w perturbation velocity component in the direction of the z axis 

x, y, z Cartesian coordinates of an arbitrary point 

xi, yi Cartesian coordinates of sovirce or doublet position 


angle of attack 


'|l-Mo" 


3m 
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Pq density in free stream 

J 

,a,ap doublet weighting factors > 

(p pertirrbation velocity potential 

Ap cpu - 

* ■ 

Subscripts 

I conditions on lower surface of wing (at z = 0—) 

u conditions on upper surface of wing ( at z = 0+) 

THE IHTEGRAL EQUATIONS 

Subsonic 


Triangular plan form .— A general solution of Laplace's equation 
which is suited to problems in linearized subsonic wing theory (given 
e.g., in reference 6) is that which relates a. velocity potential or per- 
turbation velocity to the value of its jump across a given surface. For 
a lifting triangular wing as shown in the sketch this can be written 


u = 


kit 


1 


Au dyi 


-^nxi 


[ ( x-xi ) ( y-yi ) ^+P^z^ ] 


( 1 ) 



where P \ , u is the 

perturbation velocity parallel to 
the X axis and Au is the junrp 
in u over the wing plan form. 

In linearized theory this junip car 
be related to the loading coeffi- 
cient Ap/q by the equation 



Further, the velocity potential cp 
can be found by the relation 

cp = /^ u dx 
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Operating on equation (l) in this manner and interchanging the order of 
integration gives 


cp = 



oinxx 

-illXi 


(^) 

[(y-7i)^+z^] 



X— Xi 

(x-xi )^+p2z2+32( y-7i 



which represents, physically, a distribution of elementary horseshoe vor- 
tices. 

The effect of con5)ressibility in a linearized study of lifting— 
siirface theory can only enter through the use of 3. Setting 


or 


= 1 + 


(x-xi) 

y (x-xi )^+3^(y-yi )^+p^z^ 


(4) 


it is seen that a is the only term in equation (3) which contains 3. 
This term has an interesting interpretation in the light of the study 
which has been made at sonic speeds. At Mo=l(i.e.,3=0), a has 
either the value 2 or 0, depending on whether xi is less or greater 
than X. Hence, equation (3) becomes 





^ oinxi 


Mnxi 



(y-yi)®+z^ 


(5) 


Now reversing the order of integration and using the definition implied 
by equation (2), namely, 

Ap _ _2_ 

*1 ^o 

gives finally 



A 9 dyi 

(y-yi)^z^ 


( 6 ) 


Eqxiation (6) has been studied in reference 5 as the fundamental equation 
for slender wings or wings flying at near sonic speeds. It is an equa- 
tion which gives the solution for the velocity potential in a three- 
dimensional flow in terms of two-dimensional doublets, the two dimensions 
being at right angles to the free-stream direction. A solution of such a 
nature is Immediately in^jlied by the physical character of both sonic 
wing theory, in which the Mach cone has degenerated to a Mach plane, and 
slender-wing theory, in which the wing is so slender that the chordwise 
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gradient of velocities can be neglected compared to the vertical and 
lateral gradients. 

By comparison of equation (5) with equation (3)^ ib is seen that 
the term a can be interpreted as a factor which corrects the slender— 
wing-theory results as given by equation (5) for the effects of chord- 
wise gradients in velocity and con 5 )ressibility. By consideration of the 
effect at one point of the distribution of doiiblets over the wing, this 
correction can be visualized as a reweighting of the two-dimensional 
doublets according to their position relative to the point. The sketch 
indicates the variation of a across the span at various chord stations 
for 3 = 0.6. Observe that the 
doublets ahead of the point at 
which the potential is to be 
determined are still weighted 
far more heavily than those 
behind the point. The effect 
of considering P different 
from zero, however, is to 
reduce the extreme difference 
in weight occasioned at 3=0 
so that the doxiblets behind a 
given point do have some effect 
on the induced velocities there, 
and the doublets ahead of a 
point induce a somewhat smaller 
disturbance than before. Since 
the strength of these weighted 
two-dimensional doiiblets is given 

by the magnitude of the three-dimensional loading, their' strength is 
zero everywhere off the wing plan form including the area behind the 
wing occupied by the vortex wake. 

Two different methods for the further reduction of equation (3) will 
be considered. The first method involves finding the vertical induced 
velocity for points along the x axis, while the second involves finding 
the average vertical induced velocity along the span at a given chord 
station. The first method must be discarded for triangular wings because 
of difficulties around the apex; the second, however, proves to be quite 
satisfactory. The simplification obtained by considering the vertical 
induced velocity for points along the x axis will be considered later 
in connection with the rectangular wing. 

Since it is easier to consider first the averaging process, the 
operator 
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is applied to the weighted doublets, <Jz/ [(y— yi)^+z^ ], of equation (3) 
with the result that 



where w is the average value of the verticeil induced velocity along a 
given span. 

The solution for Ap/q obtained from slender-wing theory can be 
written^ 



where in that theory fi(xi/co) = 1. If the value of Ap/q given by 
equation (8) is placed in equation (7), the resulting integral equation 
can be written in a simplified form if it is noted that 


y(x-xi)^+P^(mx-yi)^ y(x-xi)^+P^(mx+yi)2 


iiix-yi 


mx+yi 




m(x+Xi) y 

m( x-xx ) n y-m^(x^-xi^) +2Tjmx-Hi^ 


(9) 


where for the first term in the brackets the transformation q = mx— yi 
was axsed and for the second the transformation q = mx+yi. Hence, 
equation (7) finally reduces to the followin g 



The solution of equation (lO) will be discussed in a later section 
devoted to triangular wings. 

^This solution follows from an analysis of equation (6). See reference 5* 
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Rectangular plan form. — If the plan form of the wing is rectangular 
as shown in the sketch, then equation (3) is modified slightly to the 
form 

zv,^ ^ ^ 

m _ 9. / / 1.9: / 1 1 \ 


^ (y-7x)^+z" 


j-*— 2s — «-j 



It is possible in this case to 
study the vertical induced 
velocity for points along the x 
axis; that is, to find 8cp/8z 
hy equation .( 11 ) and then set 
both y and z equal to zero. 
In order to do this a specieil 
notation must be employed. Thus, 
if the indefinite integral of 
f(y)/y^ cein be written (where 
f(y) is boimded at y = O) 

f' <^7 = ^(y)+c 


Then, by definition, 


dy = J(s) -J(-s) 


By means of this definition it cein be shown that (reference 7) 






(x-xi)^+3Vi^ 


and if^ 


<1 “ ^ 0^0 


( 4 ) 


"The solution for the rectaiigular wing in slender-wing theory is that the 
load be zero across every spanwise strip aft of the leading edge. To 
find the chordwise correction to such a theory, therefore, a spanwise 
distribution must be assumed. Since, however, slender-wing theory eilso 
requires an elliptical span loading for the boundary conditions of a 
rectangular wing to be satisfied, a reasonable choice is that given by 
equation ( l4) . 
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Then, since integrating by parts gives 

T f \2 r® <^2 j 

I2 = - {x-xxrf — ^ — 

J_s 71 7 (x-; 


s^-7i^ 




equation (l 3 ) can be written 



dyi 

( 15 ) 


V 





dxx 


(16) 


This integral equation has been derived .previously by K, Wieghardt ' 
(reference 8 ) with regard to the re ctangular-^ing problem. The solution 
of equation ( 16 ) will be discussed in a later section devoted to rectem— 
gular wings. 


Supersonic 


Triangu lar plan form .— In passing from subsonic to supersonic 
theory, we pass from the. elliptic to the hyperbolic partial differential 
eqmtion and in particialar from Laplace* s eqixation to the wave equation. 
The solution which relates the perturbation velocity u at any point in 
the field to the loading on the wing can again be written in terms of an 
elementary horseshoe vortex distribution over the wing plan form. As in 
reference 9 this becomes 


-tXH 


(x-xi) Au dxi dyi 


T [ ( 7-71 ) ] J ( X-Xl ) y-yj^ ) 


and since u = ^ and — = 


9 


rr 


^ dxi dyi 


X— Xi 


(7-71 )^+Z^ 1- y(x-Xi)2-j32(y_y^)2_|322.2 


( 17 ) 


where t is the area on the wing bounded by the edges and the trace of 
the Mach forecohe from the point x,y,z. Again the effect of compressi- 
bility appears only in the term within the braces so that by defining 


Up = 


X— Xi 


J (x-xi)^-P (y-yx)^-p^z^ 


( 18 ) 
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(Jp contains all of the Mach nimiber effects at supersonic speeds. At 
Mo = 1, Op = 1, and since, by the definition of t, xi < x, it follows 
that at sonic speeds equation (17) also reduces to equation (- 5 ). Hence, 
the doublet distributions represented by equations (17) and (3) are con- 
sistent at the speed of sound. 

In order that an exact parallel can be provided with the subsonic 
solution to the triangular wing, the average vertical induced velocity 
for points along a given span is again considered. It can be shown 
(reference 7 ) that in the plane of the wing 



where the order of integration must be indicated (i.e., the integration 
with respect to yi must be made first). The letters E.P. indicate the 
real part of the term is to be taken. Such a device can be used since 
the double integral must always be a pure real qmntity in the area t 
(A p/q is real everywhere on the plan form) and a pure imaginary quantity 
over the rest of the area indicated ( see sketch) . . The average vertical 
induced velocity along the span may be obtained by applying to equation 
(19) the operator 



I3 has a derivation similar to that used for equation (9) and can be 
expressed in the final form as 
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"3 




[m(xi+x)^] [m(zi-x)+Tj] 


( 22 ) 


It is possible to find an exact solution for by means of 

equation (2l), but the discussion of this analysis is reserved for a 
subsequent section. 

Bectahgular -plein form .— Equation (19) can also be used in the case 
of a rectangular plan form by an appropriate change in limits; thus. 



where again it should be stressed that the order, of integration cannot be 
reversed. The region of integration and position of the wing with refer- 
ence to the axis is shown in the sketch. As in the case of the siabsonic 
rectangular wing, the value of w will be obtained only along the x 
axis so that y in equation (23) can be set eqiaal to zero. The loading 
will be assumed to have a form 



which is similar to that used in the 
sid) sonic case except that the refer- 
ence length is now the semi span 
instead of the chord. Such a dif- 
ference is reasonable since in the 
stq>ersonic case the position of 
the trailing edge cannot effect 
the loading on the wing. , pure 

imaginary 

Finally, therefore, when y = 0 
equation (23) becomes 




where is given by the equation 
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I = R.P. 



— s 




dyi 


(26) 


The solution to equation (25) is deferred to a subsequent section. 


LOADING ON WINGS 


The previous section was devoted to the development of the integral 
equations which are to he studied for the two types of plan forms in 
subsonic and supersonic flight. In order that this study can proceed in 
a natural manner, the arrangement of the presentation has been chafed so 
that the plan form is the principal division and the speed is subsidiary. 


Triangular Wings 


Supersonic case.- The decision to solve for the loading on the 
supersonic, triangular, flat plate by analyzing eqi^tion '(2l) was not an 
obvious one since the exact solution of the linearized partial differeit™ 
tial equation for this case has already been obtained. (See, e.g., refer- 
ences 10, 11, and 12.) Thus it is known before starting that the value 
of fs(x) in equation (2l) must be i/e where E is the con 5 )lete 
elliptic integral of the second kind with modulus However, 

these solutions were obtained by an entirely different procedure so that 
by solving equation (2l) and comparing the two results a check on the 
accuracy of the 'method xs obtained. Furthermore, when the subsonic 
problem is analyzed the same general procedure will be followed and the 
results can then be accepted with greater confidence. 

The first step in the solution of equation (2l), in which w has 
been set equal to Wg since the wing is a flat plate, is to change 
variables by the transformation = xi/x. This gives 


1 = f (x) + i igU^) 

2 « Jo 1-61 


(27) 


In the equation for 1 3 , the transformation mni = tj/x was used so that 

y (l-li)^-P" 


-L 


,22.2 

m qi 


1—^1 qi y ( 1 1 +i~q 1 ) (1 1 “1+qi ) 


dqi 
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which is con5>letely independent of x. The partial derivative of both 
sides of eq^lation (2?) with respect to x gives 


=-'■ 


jt^Pm 


l-li 


(28) 


Equation (28) is a homogeneous lineeu: integral equation. The solution 
to equation (28) is simply f3*(x) = 0 or, what is equivalent, fa(x) 
equals a constant, (^3)3 say. By means of equation (27), this constant 
can be evaluated. Hence, 


(^)o = 


nPm 




' sii.,(ei)dgi 

i-«i 


n-i 


(29) 


which represents the solution to the problem. The integral Ig was 
calcxilated analytically as in appendix A, and then the value of (fg)Q, 
as given by equation (29), was determined -by nianerical integration. 

For Pm = 0.8 the result of this computation was O.708; whereas the true 
value given by 1/E is O.705. 

Equation (27) can also be solved when the wing is slender with 
respect to the Mach cone by considering Pm to be small. Setting 
Pm = 0 yields 


^^3^Pm=0 




dll] 


ni y ( I i+i-^i ) ( 5 1-1+H 1 ) 

and this is readily evaluated to give 

, , _ 

The integral equation reduces to 

I - ^i^3(^^i)dgi 

~1 


which by a retransformation of variables xi = xg^^ becomes 

xifgCxi) dxi 


X = 


y x2— Xi2 , 


(30) 


(31) 
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Equation (3l) is a special form of Abel*s integral equation, the unique 
inversion® of which is, in this case, fg(x) = !• This is easily veri“ 
fied by direct substitution. 

The simplicity of this result is not accidental, of course, since 
the value of fgCx) was originally introduced by equation (20) as a 
correction factor to the slender-wing-theory solution. 

Subsonic case .- The study of the triangular wing presented in the 
preceding section was made first at arbitrary supersonic Mach numbers 
nntl then at a Mach ntmiber equal to 1. In keeping with this order of 
decreasing speed, the subsonic flat plate will be studied first at sonic 
speed and then for general subsonic Mach numbers. 

An inspection of equations (9) and (22) is sufficient to show that 
(li)pm^o is equal to (l3)pm=o* ^ence, equation (30) can be siibstituted 
into^equation (lO) and there resTjlts (since again w is set equal to Wq) 



and this reduces immediately to 



It is now obvious that equation (3l)^ which was derived from super- 
sonic wing theory, and eqijation (32), which was derived from subsonic 
wing theory, are identical. Clearly this establishes the continuity of 
the theoiy in. passing from the supersonic -oo the subsonic regimes. 

®If Abel’s equation is written in the form 



its inversion is 
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The analysis for the general subsonic case leads eventually to the 
numerical solution of an integral eqiiation. However, an idea of the 
qualitative form which this solution 'must assume can be gained by a 
rough preliminary analysis before resorting to the more tedious pro- 
cedures necessary for a quantitative evaluation. 


The evaluation of Ii is given in appendix A, and a plot of fcx 
against xi/x for (3m)2 equal to 0, 0.05, 0.10, and 0.20 is shown in 
figure 1. Subtracting from Ii its value at 3m = 0, and denoting this 
difference by In, it follows that 



As a rough approximation, since Iii/2n is small and tends to zero as 
3m goes to zero, consider € = Ixi/2jt to be a small constant which cein 
be. taken through the integral sign so that equation ( Slf) becomes 


X 



X Xxfl 


(il) 


dxx 


</ x2— Xx^ 



(35) 


For e equal to zero the solution to equation (35) is, of course, 
fi(xo)=l as has already been discussed. For a small but finite value 
of e, the second term will not have much' effect on fx(xo) near the 
center of the wing chord, but neeir the apex and trailing edge it becomes 
do minan t since the value fx(Xo)’ = 1 makes this term logarithmically 
infinite at these two extremes. Thus, the second term must certainly be 
reckoned with In finding the solution for fx(xo) even though e is 
small. 


Now it is apparent that the right-hand term of equation (35) makes 
that expression a singular integral equation. Experience with the singu- 
leir integr€il equation 

g(yi)dyi' 




y-yi 
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which arises in siibsonic lifting-line theory and two-dimensional airfoil- 
section theory proves useful in the present problem. ' Thus, the solution 
to equation (35) is not unique unless the Kutta condition is specified 
at the trailing edge, and, furthermore, the value of fi(xo) tends to 
infinity as x approaches the apex and to zero as x approaches the 
trailing edge. 

One would, therefore, anticipate the shape of the solution for 
fi(xo) “to something like that shown in the sketch. Such qualitative 
knowledge is useful in setting up the numerical procedure used for the 

correct solution to eqvia— 
tion (lO). The presentation 
of this procedure appears in 
appendix B and the results of 
the analysis for Aj. equal to 
0, 0.90, 1.26, and 1.79 are 
given graphically by figure 2. 

The figure shows that the 
shape of fi(xo) which was 
sketched was fairly accurate. 
Further, the values in the 
figure indicate that the center 
of pressure shifts forward with 
increasing values of ^m. A 
more comprehensive discussion 
of the integrated valiies of the 
loading will be given in a 
later section. 


Rectangular Wings 


The discussion of the triangiilar wing was divided according to the 
Mach nuBiber, The same division will be used for this section starting 
with the discussion of the resiilts for supersonic speeds, then the results 
of both supersonic and subsonic theories at so'iic speeds and f in a l ly a 
discussion of the siibsonic development. 

Supersonic case .- The solution of equation (25) will give the load- 
ing on a rectangular wing flying at a supersonic Mach number. The evalu- 
ation of the integral is carried out in appendix A where it is shown 

that I4 can be expressed in terms of complete elliptic integrals of the 
first second kinds. Having the expression for I4, a numericsQ. 

solution may be obtained for f4(x/sp) when w = Wq, that is when the 
wing is flat. (See appendix B.) Figvire 3 shows a plot of fj (a factor 
representative of the chord lift distribution) as a function of x/sp^ the 
ratio of the distemce back from the leading edge to the magnitude of the 
reduced semispan. The value of f^ given by equation (25) can be checked 
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in the interval 0<(x/sP)<2 because the exact solution to the conrplete 
linearized partial differential equation can be readily obtained there. . 
The comparison is given in figure 3. The fairly rough agreement shown 
is not surprising since equation (25) is derived on the assumption that 
the spanwise distribution of load is elliptical at every chord station, 
and certainly this assumption is least accurate in the interval where the 
comparison with the exact results is made. The area mder the exact and 
approximate curves in figxjre 3f between the initial value and that at 
which f 4 = 0, is nearly the same. (See the next section on Aerodynamic 
Characteristics.) The integrated value of f 4 as given by equation (25) 
therefore, should be sufficiently accurate for x/s^> 2. 

As for the qualitative nature of the variation, figure 3 shows that 
the loading on a narrow rectangular wing flying at supersonic speeds 
falls linearly to zero, becomes negative, and then oscillates between 
negative and positive values; the amplitude of the oscillation being so 
heavily damped that after the third change in sign the magnitude is 
practically zero. 


It should be noticed in studying the results of figure 3 that the 
entire resultant lift of the wing is concentrated in the interval ' 
0<(x/sP)<2. But as the Mach n\miber approaches 1 this interval 
approaches zero, and the entire lift of the wing is carried in a strip 
along the leading edge.^ Such a solution violates, in the vicinity, of 
the leading edge, the assumption on which the theory is based an d should 
be considered only as a theoretical limit. ^ 

Be suits for the lift and pitching moment on the recteingular wing 
will be developed in a later section. . 

Subsonic case .— The study of the subsonic rectanguletr wing stems 
from equation (l6). The first step in the analysis of the equation will 
be to consider its solution at 3s = 0 and show that this is continuous 
with the supersonic re stilts there. 


The value of Is can be written (equation (l5)) as 



and for 3s = 0 this becomes 


'1 /Cx-xi ) ^+3^ s^y^ ^ 

= / / 

(x— xi)n xi <x 

— (x— xi)rt xi >x 


Ip = 


'^his result also follows by inspecting equation (25) for the values 
3s = 0. 


1 ' 
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and hence equation (l6) can be written 

( 36 ) 

Equation. (36) is identical with the form of the supersonic equation (25) 
at Ps = 0 so that once again the continuity of the subsonic and super- 
sonic theories at the sonic speed range is established. Furthermore, 
equation (36) shows that if w/wq is constant then f2(xo) must be zero 
everywhere except at points where it can be represented by a pulse' the 
integral of which has a finite magnitude. From the supersonic discussion, 
it is clear that one such pulse exists and is located at the leading edge. 

The evaluation of I2 for 3s >0 is given in appendix A. The 
nimierical solution to equation (16), assiiming the Kutta condition at the 
trailing edge, is given in appendix B for values of reduced aspect ratio 
Ap (defined as 3 times aspect ratio) equal to 0.33^ 1.0) 1.5> and 2.0. 
For an aspect ratio equal to 2, these values correspond to Mach mmibers 
of 0.986, 0.866, 0.662, and 0,. respectively. The results of the compu— ^ 
tations are shown in figure 4 where the chordwise lift distribution 
factor f2(xo) is plotted against Xq for the various values of Aj.. 

By comparison of figure 4 with figure 3> it can be seen that in the 
subsonic case the loading drops monotonically from infinity at the lead- 
ing edge to zero at the trailing edge and does not oscillate in the after 
portion, as in the case of the supersonic wing. 

When 3 equals one, these results can be compared with those 
obtained by Wieghardt and presented in reference 8. The sketch shows 
the comparison for two values of the aspect ratio. Curves are also 
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shown in the figure for the loading obtained by using the method given 
in the appendix of this report but by satisfying the integral equation 
at only six and three points. The latter curve is in better agreement 
with Wieghardt's result and, since Wieghardt (although using a differ- 
ent method involving Birnbaum functions) used only 4 points, this may 
accoxmt for the discrepancy between the final results of this report and 
those of Wieghardt. 


AE310DTNAMIC CHABACTEEISTICS 


The previous section presented solutions for the loading on trian- 
gular and rectangular wings flying at subsonic and supersonic speeds. 
This section will be devoted to the conversion of these loadings to 
expressions for lift and center of pressure. 


Lift 


By definition the lift coefficient can be written 

= 1/^ f dy (37) 


and this will be evaluated for the various cases for which the loading 
coefficient has been obtained. 


Supersonic triangular wing .— Since the exact linearized value for 
the loading on the triangular wing flying at supersonic speeds has been 
derived, the lift coefficient can be written in the form 



(38) 


where A is the aspect ratio and E 


second kind with modulus 


k 


is the elliptic integral 'of the 


S\ 3 b sonic triangular wing .— In the case of the subsonic triangular 
wing, equation (37) becomes 


,Co 


mco^ Jo . 

and this becomes (since A = 4m) 


^m^x 


— mx 
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Ar ’.90- 


oA 


Af^i79 


fi(xo)^c 


The numerical evaluation of equa- 
tion (39) is not difficiilt since 
xofi(xo) has the variation 
shown in the sketch. 


Supersonic rectangular wing ,— For values of Cq<2sP the exact 
value of the lift coefficient on a rectangular wing flying at supersonic 
speeds has been obtained and can be written in the form: 


For Aj. > 1 


Cl 

oA A-r^ 


When Ap<l equation (3?) must be used in connection with equation (2k) 
and there resiilts 


Cl = r ° dx r — f 4 Vs^-y^ dy 

2scoJo V-s s3 \sp/^ 


which reduces to 




this can be written in the form 


oA Aj-Cq Jq V ' 


which becomes, if xa = — —y 

CoAr 


for A-r < 1 


Cl ^ 

oA 2 Jo 


f^ (xa) dxa 


SubRonic ar wtng .- The equation for the loading on a sub- 

sonic rectangular wing, equation (l^)> placed in the formula for lift 
coefficient yields 
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which becomes 


9l = 

oA 



(42) 


Th© evalua’tion of equation (4-2) ty numerical means requires special 
consideration since faC^o) approaches infinity at the leading edge as 
shown in figure 4. To this end, rewrite equation (42) in the form 


^ " 2 X ^ 2 (^ 0 ) <^o ^ 2 (43) 

and equation (b 4) in the appendix (for the special case in which x^=l) 
in the form ° 


1 = ^ r f2(xo) 

2jt J Q 


It + 


At® 
k(l-xo) J 


dxo + 


( 44 ) 


An application of the mean— value theorem yields 



n+— 

-k0(i-0)j 



Tt + 



( 45 ) 


where Eg has the modulus ke which equals, Aj./y4(l-02')+,Ai.2 and where 
0< 0<€. The combination of equations (43) and (45) yields an expression 
for the lift coefficient involving only the load distribution from a 
distance e/co back of the leading edge to the trailing edge. 


Pitching Moment 


By definition the pitchlng-^noment coefficient about the apex or 
leading edge and based on the root chord can be written 
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Cm = - ^ f f X ^ dxdy (^6) 

Sco J Js <1 

Equation (i|-6) will be applied to the various loadings which have been 
studied. 

Supersonic triangular wing .- The exact linearized value for the 
pitching-^noment coefficient on a triangilLar wing flying at supersonic 
speeds has been derived elsewhere and can be written in the form 

^ = — 5- (l<-7) 

oA 3E 

Subsonic triftpgular wing .— The derivation of the pitching moment 
on a subsonic triangular wing proceeds in the same manne r as the deriva- 
tion of lift emd there resiilts 


Hffl = - jt r xo^fi(xo) dxo 

1 


m 


This expression can be easily integrated numerically since XQ^fi(X(. 


varies as the sketch indicates. 



Supersonic rectangular wing .— 
For values of reduced aspect ratio 
A^ greater than 1 the pitching- 
moment coefficient on a rectangular 
wing is given by the equation. 


for Ar> 1 


oA 


6Aj» '^1 

3Ar2 


(49) 


When Aj.< 1 the solution to the 
integral equation must be used and 
the final expression can be written 


a/Ar 


xaf4 (xa) dxa 


(50) 


Subsonic rectangular wing .— The equation for the pitching-nmoment 
coefficient on a subsonic rectangular wing follows in the same manner as 
did that for the lift coefficient. 
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Hence, 

- = Xof 2 (xo)dxo (51) 

oA 2 Jo 

and, since the variation of 
XofaCxo) is as indicated in 
the sketch, the numerical 
integration of equation (5l) 
is simple. 


Center 



of Pressiare 


Since the pitching moment is based on the root chord, the center of 
press\ire of all wing plan forms can be written 



.p. 


o 



(52) 


Discussion of Eesults 


Figures 5 and 6 show the variation of the lift coefficient and 
center of press\xre for triangular and rectangular wings for the range in 
which the reduced aspect ratio Aj. is small. For the triangular, wing 
the differences between the subsonic and supersonic cases are not large 
in the interval of A^ shown; the siibsonic wing develops somewhat less 
lift and the center of pressure moves forward as Aj. increases. The 
characteristics of the rectangular wing, however, show a large variation 
in passing through the speed of sound. 

The subsonic rectangular wing has a variation of Cj/Aa with Ay 
which is almost identical with that for the subsonic triangulsir wing. 
.Unlike the triangular wing, however, the curve for p /c© on the 
rectangular wing shows that this lift is carried farther and farther 
forward from the quarter-chord position at Mp= 0 all the way to the 
leading edge at M© = 1. 

As the speed is further increased and the rectangular wing enters 
the supersonic speed range, the magnitude of the lift oscillates about 
the curve for the subsonic case until a reduced aspect ratio of about 
O ..5 is reached and then rises to a maximum at about the point where Aj. 
equals 1. The exact curve (obtained by linearized lifting-surface ■ 
theory) for values of Ap greater t’nan one is shown as a dotted line for 
the purpose of comparison. 
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The variation of the center of pressure on a supersonic re ctangul^ 
wing indicates that the wing is unstable for. all 
pSSt behind the leading edge for values of Ar 

of pressure, in such a range, having moved forward of the wing leading 
As Ar increases past the value 0.5^ however, the center of 
pressure moves back along the wing and rapidly approaches the midchor 
point, its location according to two^imensional linearized theory. 


C0NCLI3DING REMAREB 


A theoretical investigation has been made of the correction to 
slender-^lpg theory which accounts for the effects of compressibility 
and the chordwise gradients in velocity. Integral equations were devel 
oped and numerically evaluated for triangular and rectangular plan forms 
in both subsonic and supersonic flow. 

In the case of the triangular plan form, the value of the lift coef- 
ficient predicted by slender-wing theory was reduced by ap which 

increased with increasing aspect ratio as was to be 

reduced aspect ratio of 1.0, this correction amounted to 6 percent for 
both subsonic and supersonic flow, but, at a reduced aspect ratio of 
Lo, the correction had increased to 21 percent in the case of subsonic 
flow and 1? percent in the case of supersonic flow. 

The center of pressure of a triangular wing in supersonic flow 
remained constant as predicted by slender-wing theory; in subsonic flow, 
however, there was a slight forward shifting of the center of pressure 
with increasing aspect ratio which amounted to 6 percent of a root chord 
length at a reduced aspect ratio of 2.0. 

In the case of the rectangular plan form in subsonic 
rection to the lift coefficient due to a consideration of the chordwise 
and compressibility effects was the same as for the triangilar plan form. 
In supersonic flow, below a reduced aspect ratio of 0.5, thj 
lift coefficient oscillated about the constant value predicted by slender 
wing theory. Above this, value the lift coefficient was 
ing a maximum correction of 32 percent at a reduced aspect ratio of 1.0. 
At Aj. = 2.0 the lift coefficient was again decreased by 9 percent. 

According to slender-wing theory, the center of pressure of a 
tangular plan form is located on the leading edge. In the subsonic case, 
tS^rrections shifted the position of the center 

the leading edge as the aspect ratio was increased. This shift amomted 
to 19 percent of a chord length at a reduced aspect ratio of 2.0. in 
the Lpersonic case, a consideration of the chordwise and compressibility 
effects indicated that the wing was unstable for a pivot point located 
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the leading edge around a reduced aspect ratio of 0»4, As the aspect 
ratio increased, the center of pressinre moved back of the leading edge 
toward the midchord position. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Nov. 28, 1950. 
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APPENDIX A 

EVALUATION OF SPECIAL . INTEGRALS 
THE INTEGRAL Ii 


The evaluation of Ii will he discussed first for the case in which 
x>xi and second for the case in which x<xi. 

Case 1, x>xi 


It is possible to write 1^ in the form 

pP-i J 

Ii = m / ■ ■ o-n 


(Al) 


where |io = m(x-xi) and [xx = m(x+xi). Th,e linear term in the lower 
radical of the integrand can be eliminated by the transformation 
= (d + 6t)/(l + tT. and the integral becomes 


Ii = 


m(5-<T) - g-^i y 


/ M’l-S r— r 

vO^o-^Ko^ ( i+t) (0+6t) A- ( t‘ 


dt (A2) 


where 


a = 


~^o( ) -H^o y ( 

P^m^(^i+|ao) 


(A3) 


and 


& = 


y( 1+P^m^ ) ( ) 


(Ai+) 


The expression for a and 6 may be combined to give the useful identi- 
ties 


= — obp^m^ 

(6-Pq) = 0 
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Using fundamental properties of even and odd functions, equation (A2) 
may be reduced to the form ' 


Ii = 2m 


6) Jq 


£ 


S/q 


1+ ur 1+ cj 

6k ak*^ 


k« 


'k «^+k^fa)^ 
1-0)2 


do) 


by the substitution 


0 ) = 


hi-S 

Pi-<j 


and where 


2 (m-<T)^8 


k‘=' = 


(pi2-CT8)(8-cr) 


k2 


_ 8 / 
k‘2 0 \ M.1-8/ 


By introducing the Jacobian elliptic functions in the transforma- 
tions 0 ) = cnu, the integral reduces to 


Ii = 


where 


^Uq2+62P2jji2 0 

Y V,- 

2. V 2 
6 2- \ 

(k^l-<j)(cT-do) "^0 


l+Va^sn^u ' 


-(§> > 

0 ^ • 0 

0 


k« + 1 k 



) dn^udu (A5) 


(A6) 
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and 




\ 

where 

G(x) =EF»(x) + KE*(x) -^KFKx) (A9) 

the modulus of the elliptic integrals being k or k». 
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a = 


^ao(Ho-^2ln^^il)-^Xo y(l+^2In^)(^Xo2+^2m^^ll2) 




(All) 


In this case it can be shown that in equation (A5) V^> 0 and 


V < —1, and the solution for Ii is 
2 


Ii = 2m 


{K T^y (^)^(^) ] 


1 V / 

8 y (|iiHJ)(cr-do) 



where G(x) is defined as in equation (A9), 


THE HJTEGEIAL I2 


Writing I2 in bbe form 


and setting 



I2 can be integrated to give 


I2 



2PsE 

k 


THE IMEGEAL 


(A12) 


(AI3) 


(Al4) 


Case 1, 0o^^i< iio< 
Writing I3 in the form 
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I. = e. 






and making "bliB transfonuation 




where 


2Qq ^tl - 
1— Oq ^^o 


reduces equation (A15) to the form 


/ 2^lQ / oK dll ^pdii 

0o(^ll-^lo) \^o 1+ k^sn^u 1- k2sn2u 


The integration may now he completed and 


^lx•Hlo 


k \ (l— + 



^l-0n^ G 


Vl + k^ 


2 




Case 2, 0 < Po '^0o^ll 


In this case 


^o ^ 
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The integration can be completed so that 




THE INTEGRAL 
Case 1, 0^i^x-3s 
The integral- can be written 



(A20) 


(A21) 
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When, the transformation snu = yi/s is made, the expression may be 
integrated to give 

pK 

14 = 2 (x— xi) / dn^udu = 2{x— X i)Ei (A22) 


where 


ki2 = 


2o2 




(x-xi)' 


(A23) 


Case 2, x-Ps<xi<x 
In this case 1 4 can be written 




(x-Xj^)^-P^yj^ 
s2-yi2 




(A2k) 


The transformation snu = Pyi/(x— xi) applied to eqioation (A24) yields 

>K (x— Xj^)^ 


1a = 2 


Ps 


cn^udu = 2 k 2 (x— xi)B 2 


where 


and 


k2. = 


X— X, 


Ps 




(A25) 


(A26) 


(A27) 
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APEENDIX B 

NUMEEICAL SOLUTION OF INTEGRAL EQUATIONS 
SUBSONIC TRIANGULAR WING 


Since the integral Ii is a function only of the ratio xi/x, 
equation (lO) can he written for w = Wq 



where | = xi/x. It is now assumed that f(x|/cQ) may be considered 
consteint over smaJl intervals. This reduces the solution of the integral 
equation to the elementary problem of solving a system of simultaneous 
algebraic eqmtions. On the basis of such an assumption equation (Bl) 
becomes 


-27 {I 


where 


3 


go(|) = -pSi=+ — li(0 0<| <1 ' 

Vl-r 1-1 

60(5) = A K|<V 

X 5 X / 


(B2) 

• • n 


(B3) 


The function go(0 can be integrated by numerical means, and it is 
convenient to make a plot such as is shown in the sketch in which each 
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ordinate represents the value of go(0 integrated over the interval 


- .90 


Systems of simultaneous equations were obtained for values of n 
equal to 3^ .6, and 9» Solutions were found using the Gauss-Seidel 
method (for which the simultaneous equations were well suited) and the 
convergence was rapid as the sketch indicates. 


o J points 
A 6 points 
a 9 points 


0 


.2 


4 


.6 
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SUBSONIC RECTANGULAR WING 


Sxibstituting the value of I 2 given by equation (Alif) into equa- 
tion (16), one has for w = Wq 


1 = — 
- 2it 


f 2 d ^ 


where 


and where 


x-xi X _ AyE 

^r\ ) ^ v—'v 1 


TO 


(x-xi)2+p2g2 y4/^y + 


-A satisfactory numerical solution of equation (b 4) requires the solution, 
of the system of simtiltaneous equations of the form 




„ f2\^ /2i-l 2\ 

• ^(ny®( 2n nj"^ 


-¥) ] 


fa(l)g - 1 
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The convergence of the solutions to equation (B5) is indicated in the 
sketch where the value of n was successively taken to* be 3^ 6, and 9. 


/.2 


f, 1.0 


— 


AO- 


•OA— A 


o 3 points 
A 6 points 
a 9 points 




□A O 


'x A 


A -1.26 


SUPERSONIC RECTANGULAR WING 


Equation (25) can be written when w = w^ as 


1 = f. 




(t 


where from equations (A22) and (A25) 


^>1 

Jt Ps 


2 Sii < 1 

tT ^ ^ ft c — 


and where 
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X— Xi 


X— Xi 

sP 


By the application of the trapezoidal rule for nvunerical integration to 
equation (b 6), it, is possible to write f^Cx/sP) explicitly as 



where A(x/Ps) is the inte^al of the trapezoid. 


(B8) 
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Figure /.- Variation of j, with x,/x . 
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Figure 2.~ Variation of chordwise correction ^ factor f^ 
for subsonic triangular wing. x^fffx^J dxc, 
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Figure 4,— Variation 
f£ for Bubsonic 


of chordvtiBo eorroction 
rectangular wing. -JiAf'fpfxol 


factor 
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(a) Lift. 



Figure 5.- Aerodynamic clioracterietics of a triongular 
wing having a low value of Am . 
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{q) Lift, 


^c.p. 

Co 



(b) Center of Pressure, 


Figure 6.- Aerodynamic characteristics of rectangular 
wing having a low value of Af-. 
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